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Several  continuing  investigations  of  the  seismic  source  us.ng 
near-  and  far-field  data  are  discussed.  We  completed  our  analysis  of  near- 
field acceleration  data  recorded  at  Bear  Valley,  California,  by  examining 
two  small  events  which  occurred  on  7 February  1974  and  6 Tuly  1974.  The 
February  event,  using  only  spectral  analysis,  was  found  to  have  a moment 
of  3 x 1021  dyne-cm  with  an  equivalent  circular  dislocation  ladius  of  0.  18 
km2.  Both  spectral  and  time  domain  waveform  fitting  techniques  were  ap- 
plied to  the  July  event.  The  solutions  using  each  method  were  in  reasonable 
agreement,  with  a seismic  moment  in  the  range  of  2 to  4 x 10  dyn,-cm 
over  an  equivalent  dislocation  radius  of  about  0.  1 to  0.  3 km  . 

Using  far-field  surface  wave  data,  two  central  California  area 
earthquakes  (22  June  1973  and  28  November  1974)  were  re-examined  using 
two  new  earth  models  which  were  determined  by  McEvilly  (1975)  for  the  north 
east  and  southwest  side  of  the  San  Andreas  fault.  Several  far-field  solutions 
for  the  source  mechanism  were  obtained  by  using  various  combinations  of 
stations  and  earth  models,  including  use  of  both  structures  at  once  for  the 
first  time.  However,  for  the  central  California  event  of  ?8  November  1974, 
the  mechanism  which  agreed  closely  with  that  obtained  by  other  means  used 
only  those  stations  on  the  west  side  of  the  fault,  with  the  southwest  structure 

used  in  the  spectral  fit. 

Theoretical  first  higher  mode  surface  wave  spectra  were 
generated  for  a double  couple  source  in  a layerec  half  space.  Both  Gutenberg- 
Bullen  and  Hamilton-Healy  earth  models  were  used.  As  with  fundamental 
mode  spectra,  'holes'  were  found  almost  exclusively  for  vertical  strike-slip 
faults.  Higher  mode  spectral  levels  in  the  3 to  9 second  period  range  were 
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found  to  be  comparable  with  15  to  50  second  period  fundamental  mode  spec- 
tra. Also,  from  the  excitation  of  the  earth  models,  greater  amplitude  higher 
mode  signals  from  shallow  events  should  occur  from  high  velocity  crustal 
source  regions. 

Preliminary  results  have  been  obtained  for  five  large  events 
which  were  a part  of  the  Sinkiang,  China  earthquake  swarm  of  August  1974. 
Both  simple  spectral  and  spectral  ratio  fitting  methods  were  employed.  The 
far-field  solutions  showed  a consistent  trend  for  these  five  events  of  shallow 
focal  depth. 

Finally,  we  concluded  our  examination  of  the  scatter  of  the 
earthquake  population  in  the  M - mb  discriminant.  Focusing  on  the  body- 
wave  magnitude  estimate  instead  of  the  surface  wave  magnitude  as  in  pre- 
vious studies,  we  applied  the  maximum  likelihood  method  of  Ringdal  (1975) 
in  an  attempt  to  improve  the  accuracy  of  the  PDE  bodywave  magnitudes.  Al- 
though the  scatter  of  the  population  did  not  change  significantly,  the  slope  of 
the  data  set  decreased  approximately  0.25. 
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SECTION  1 
INTRODUCTION 

la  '.his  report,  we  present  a discussion  of  several  continuing 
investigations  of  the  seismic  source  from  near-  and  far-field  data.  We 
complete  our  analysis  of  near-field  acceleration  data  recorded  at  Bear  Valley 
in  Section  II.  Two  rather  small  events  (mb's  of  3.  0 and  3.  2)  that  were  re- 
corded in  1974  (7  February  and  6 July)  were  studied  using  both  spectral 
methods  and  time  domain  waveform  techniques. 

Several  studies  using  far-field  surface  wave  data  to  determine 
source  parameters  are  discussed  in  Section  III.  Two  central  California  area 
earthquakes,  22  June  1973  and  28  November  1974,  are  re-examined  using 
two  new  earth  models  by  McEviily  (1975)  for  the  northeast  and  southwest 
sides  of  the  San  Andreas  fault.  In  particular,  we  examine  the  problem  of 
fitting  far-field  surface  wave  spectra  when  the  source  region  contains  a 
major  discontinuity.  In  this  section,  we  continue  our  examination  of  theo- 
retical first  higher  mode  spectra  for  characteristics  that  would  reflect 
source  depth  and  structural  dependence,  such  as  'holes'  and  amplitude  level 
over  a particular  frequency  range.  Also  in  this  section,  we  present  a pre- 
liminary analysis  of  several  events  that  occurred  during  an  earthquake  swarm 
in  Sinkiang,  China,  in  August  1974. 

Finally,  in  Section  IV,  we  conclude  our  examination  of  the 

scatter  of  the  earthquake  population  in  the  M - m.  discriminant.  Instead  of 

s b 

focusing  on  the  determination  as  in  our  previous  studies,  we  apply  the 
maximum  likelihood  method  of  Ringdal  (1975)  in  an  attempt  to  improve  the 
accuracy  of  the  PDE  bouywave  magnitudes  for  our  data  base. 
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SECTION  II 

ANALYSIS  OF  TWO  SMALL  BEAJ1  VALLEY  EVENTS 
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A.  INTRODUCTION 

Since  the  occurrence  of  a magnitude  4.2  earthquake  in  June  of 
1973  two  additional  events  having  potential  for  near-zone  analysis  have  been 
recorded  by  stations  of  the  University  of  California  accelerometer  array  lo- 
cated along  the  Bear  Valley  segment  of  the  San  Andreas  Fault.  These  events 
were  a magnitude  3.  0 earthquake  which  occurred  on  February  7,  1974  and 
one  of  magnitude  3.  2 which  was  recorded  on  July  6,  1974.  While  being  of 
significantly  lower  magnitudes  than  earthquakes  previously  studied  in  the  near- 
field, accelerometer  records  of  relatively  high  signal-to-noise  ratios,  pro- 
bably accounted  for  by  the  shallow  depths  associated  with  these  events,  were 
obtained  at  a sufficient  number  of  stations  to  warrant  analysis  to  determine 
the  source  characteristics.  A.s  with  the  earlier  June  earthquake,  (Turnbull 
and  Battis,  1974  ) , the  intermediate -zone  coverage  of  the  July  6 earthquake, 
resulting  in  determinations  of  location  and  fault  plane  solution  independent  of 
near-field  observations,  along  with  the  relatively  wide  azimuthal  coverage  in 
the  near-field  would  provide  a set  of  restrictions  on  any  solution  resulting 
from  the  Haskell  moving  dislocation  source  model.  For  this  reason  and  as  a 
test  of  the  applicability  of  the  Haskell  earthquake  source  approximation  to 
events  of  very  low  magnitude,  the  near-zone  analysis  concentrated  on  the  July 
6,  1974  earthquake. 

The  University  of  California  Bear  Valley  accelerometer  array 
is  located  approximately  30  km  southeast  of  Hollister,  California  along  a seis- 
mically  active  segment  of  the  San  Andreas  Fault.  In  the  vicinity  of  the  array 
the  main  fault  zone  strikes  to  the  northwest,  separating  two  distinct  crustal 
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structures.  To  the  southwest  o(  the  fault  lies  a region  of  high  velocity  g 
a„d  metamo rphic  strata  with  tower  vetoed  sedimentary  rocks  predominating 
the  northeastern  bloc*.  Detailed  crustal  studies  have  been  conducted  for 
this  region  (McEvilly.  pe rsonal  communication.  1974;  AM  and  Lee,  1974,  ano 
Engdahl  and  Lee,  1974)  which  have  indicated  a structure  of  complex  mhomo 
gene ities  within  the  fault  rone.  As  the  Haskell  moving  dislocation  source 
model  is  a whole  space  representation,  the  velocity  profile  determined  by  Me 
Evilly,  as  given  in  Table  11-1.  was  used  to  evaluate  two  sets  ol  'whole-space 
velocities  for  use  with  stations  on  each  side  of  the  fault.  These  velocities 
were  obtained  by  using  a weighted  average  of  the  velocity  profile  from  the 
hypocentral  depth  to  the  surface  with  the  weighting  factors  for  each  velocity 
being  proportional  to  the  thickness  of  its  respective  layer.  These  weighting 
factors  and  the  whole-space  velocities  for  each  crust  block  are  also  given  ,n 
Table  11-1.  It  should  be  noted,  however,  that  in  the  case  of  spectral  evalua- 
te o,  the  source  characteristics  using  Bruno's  model  ,1970).  it  was  felt  that 

■ ■ (i,„  oarthaiiake  hvpocenter  under  study 

the  velocities  of  the  layer  containing  the  earthquake  nyp 

were  more  appropriate  for  use  in  analysis. 

In  Figure  11-1.  the  location  of  each  station  in  the  near-field  ar- 
ray, as  operational  during  the  recording  of  the  two  events  studied  in  this  re- 
port is  shown.  Each  site  is  clipped  with  a three  component  accelerometer. 
The  dashed  line  in  this  figure  represents  the  surface  trend  of  the  San  n reas 
Fault  through  Bear  Valley.  In  an  effort  to  improve  the  ammuthal  coverage  o 
events  occurring  Cose  to  the  near-field  array,  several  modifications 
network  configuration  have  been  made  since  the  June  *3.  1973 
primary  importance  is  the  shift  o,  station  6 from  the  northeast  to  the  south- 

r , Thin  added  a second  station  located  on  the 

western  side  of  the  fault  zone.  This  added  a se 

. t . au  e-mnihwpst  which  aids  in  evaluation  ol 
higher  velocity  crustal  material  to  the 

t . • _ / r pee ive r location  relative 

differences  in  signal  character  as  a func  1 

to  the  fault. 
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FIGURE  II- 1 

MAP  OF  BEAR  VALLEY  ACCELEROMETER  ARRAY  SHOWING 
EPICENTERS  OF  THE  FEBRUARY  7,  1974  (E-l)  AND 
JULY  6,  1974  (E-2)  EARTHQUAKES 
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For  this  report  the  source  parameters  of  both  the  February 
and  the  July  earthquakes  were  determined  using  spectral  techniques.  In  addi- 
tion, for  reasons  stated  earlier,  time  domain  waveform  fitting  was  used  with 
the  July  event.  In  the  following  sections  the  methods  used  in  the  analysis  and 
the  results  of  the  analysis  cf  each  event  is  presented  and  the  implications  of 
these  results  are  discussed  in  relation  to  nr.  earlier  examination  of  the  June 
22,  1973  Bear  Valley  event  (Turnbull  and  Battis,  1974). 

B.  EVALUATION  OF  SOURCE  PARAMETERS 

The  techniques  used  to  evaluate  source  parameters  from  near- 
zone data  have  been  covered  in  detail  by  earlier  reports  (Turnbull  and  Battis, 
1974;  Turnbull,  et  al.  , 1 97 S) . However,  to  facilitate  the  understanding  of 
the  results  presented  here,  a brief  discussion  of  the  two  methods  used  in  this 
report  are  given  below. 

Using  Brune's  source  model  (1970)  it  is  possible  to  estimate 
the  seismic  moment,  , and  a measure  of  the  area  of  faulting,  r , the  equi- 
valent circular  dislocation  radius,  from  the  displacement  amplitude  spectrum 
of  a given  record  of  an  earthquake.  The  equations  giving  these  relationships 
are: 


M - 47tp/33R  fl 
o o 


(II-.1) 


and 


0/<2fffc) 


(II-2) 


where  P is  the  material  density,  P the  propagation  velocity  of  the  type  of 
bodywave  dominating  the  time  gate  used  to  calculate  the  spectra,  R is  the 
hypocenter  to  station  distance,  all  supposedly  known  parameters,  and 
is  the  low  frequency  spectral  level  and  f is  the  corner  frequency  of  the 
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spectra.  1 addition,  other  parameters  of  the  source,  such  as  stress  drop, 
can  be  evaluated  from  this  model. 

In  this  study  the  acceleration  amplitude  spectrum  of  each  com- 
ponent is  calculated  directly  from  the  acceleration  data  as  a preliminary  step 
to  evaluation  of  the  associated  displacement  spectrum.  Due  to  the  superior 
noise  suppression  qualities  of  the  maximum  entropy  spectral  estimator  over 
the  discrece  fast  fourier  transform  (Turnbull,  et  al.  , 1975)  the  maximum 
entropy  method  is  used  for  this  step.  The  acceleration  amplitude  spectrum  is 
then  converted  to  a displacement  spectrum  by  two  subsequent  divisions  of  the 
spectra  by  (2  TT f).  At  this  point  all  necessary  informal  on  for  the  evaluation 
of  the  corner  frequency  and  low  frequency  level  for  a given  component  is  avail- 
able. 

The  first  division  of  the  acceleration  amplitude  spectrum  yields 
the  ve'ocity  amplitude  spectrum.  The  corner  frequency  of  the  displace  .lent 
spectrum  is  the  frequency  of  the  peak  value  of  the  velocity  amplitude  spectrum. 
As  noise  in  the  acceleration  spectrum  will  usually  cause  the  extreme  low  fre- 
quency velocity  amplitudes  to  blow  up,  these  frequencies  are  neglected  in 

evaluation  of  f . Given  the  value  of  f_  the  low  frequency  level  is  then  cal- 
c ^ 

culated  as  the  avenge  displacement  spectral  amplitude  between  the  point  at 

-2 

which  the  low  frequency  spectrum  deviates  from  a slope  of  (2 7T f ) and  the 

corner  frequency. 

To  obtain  the  values  used  for  evaluation  of  source  parameters 

the  information  from  the  components  of  each  site  must  be  combined.  In  the 

case  of  fl  this  is  done  by  resolving  the  individual  component  low  frequency 
o 

levels.  In  the  case  of  corner  frequency  the  intermediate  value  is  selected. 

This  would  seem  to  be  more  rational  than  simple  averaging  as  the  latter  may 
result  in  a v ilue  of  fc  which  could  not  be  supported  by  the  spectra.  As  was 
stated  earlier,  these  values  are  hen  used  in  the  equations  of  Brune's  model 
to  obtain  estimates  of  the  source  parameters. 
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The  second  method  of  source  charactr  ristic  evaluation  used  in 
this  report  is  Haskell's  moving  dislocation  model  (Haskell,  1969).  In  this 
process  the  time  domain  acceleration  tiaccs  are  numerically  integrated  to 
yield  velocity  waveforms  for  each  component  of  a station.  As  the  Haskell 
model  is  a whole  space  representation,  the  complex  structure  of  the  P-  and 
S-waves  packets  can  not  be  simulated  by  this  model.  Thus,  by  some  means 
a single  period  of  the  P-  or  S-wave  must  be  isolated  for  each  station  to  be 
used  in  the  fitting  process. 

Using  external  information  such  as  hypocenti  r locations  and 
fault  plane  solutions  to  provide  an  initial  case,  theoretical  waveforms  for  unit 
transverse  and  longitudinal  dislocations  of  the  fault  plane  are  calculated  for 
each  station.  At  this  point  adjustments  to  fault  length  and  width,  rise  time, 
and  rupture  velocity  are  made  to  produce  theoretical  pulse  durations  approxi- 
mating those  of  the  empirical  waveforms.  Again,  restrictions  such  as  spec- 
tral estimates  of  fault  area  ar  3 used  to  constrain  the  parameters.  Following 
these  adjustments  least  squared  error  combination  of  the  two  dislocations 
are  calculated  for  fitting  each  station  separately  and  as  a group.  The  fault 
plane  parameters  of  dip  and  strike  are  then  perturbed  to  locate  the  minimum 
residual  solution  for  each  station  and  for  all  stations  together.  Once  this  set 
of  solutions  is  determined  the  theoretical  waveforms  for  each  solution  at 
every  station  is  calcula'ed.  At  this  point  a qualitative  judgement  on  the  'good- 
ness of  fit1  is  made  to  determine  which  solution  is  acceptec 

It  should  be  noted  that  the  least  squares  error  combination  for 
all  stations  taken  as  a whole  may  not  necessarily  be  the  best  solution.  This 
situation  can  result  from  such  things  as  uncertainties  in  determining  which 
pulses  to  fit  for  some  of  the  stations  or  interference  patterns  effecting  the 
empirical  amplitudes.  Under  these  circumstances  the  observer  may  select 
one  of  the  solutions  based  on  a single  station  as  being  a superior  fit  ‘o  the 
empirical  waveforms. 
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THE  FEBRUARY  7,  1 V74  EARTHQUAKE 


On  February  7,  1974  an  earthquake  of  local  magnitude  3.0  oc- 
curred along  the  San  Andreas  Fault  in  Bear  Valley,  California  and  was  re- 
corded by  three  s’ations  of  the  University  of  California  accelerometer  array. 
The  epicenter  of  this  event  was  located  approximately  1.  3 kms  southeast  of 
Station  4 of  the  array  (McEvilly  and  Johnson,  1974).  Ihe  location  of  this  epi- 
center is  shown  in  Figure  11-1  by  the  closed  circle  labeled  E-l.  The  depth  of 
the  event  was  calculated  to  be  5.28  kms.  Considering  possible  location  errors 
and  variations  in  the  fault  surface  with  depth  it  would  appear  that  this  event 
occurred  on  the  main  fault  surface  of  the  San  Andreas  Fault. 

Probably  due  to  the  low  magnitude  of  the  earthquake,  only  the 
records  of  the  three  stations  closest  to  the  epicenter,  Sites  2,  4,  and  8,  were 
of  satisfactory  quality  for  analysis.  As  an  example  of  the  records  produced 
by  this  event  the  acceleration  and  derived  velocity  traces  from  Station  2 are 
displayed  in  Figures  11-2  and  II-3.  It  is  apparent  that  the  S-wave  portion  of 
the  signal  is  dominated  by  a much  higher  frequency  signal  than  n previously 
studied  earthquakes  (Turnbull  and  Battis,  1973;  Turnbull  and  Battis,  1974). 
However,  while  the  S-wave  packet  of  the  velocity  waveforms  examined  here- 
tofore tended  to  be  dominated  by  a single  pulse,  it  can  be  seen  in  Figure  11-3 
this  is  not  true  for  the  February  earthquake.  As  the  isolation  of  a dominant 
pulse  is  essential  to  the  successful  application  of  the  Haskell  model  it  was 
decided  to  use  only  spectral  methods  to  ascertain  the  earthquake  source  para- 
meters of  this  event. 

The  acceleration,  velocity  and  displacement  normalized  energy 
density  spectra  of  the  S-wave  estimated  by  the  maximum  entropy  method  for 
the  N45E  component  of  Station  2 are  shown  in  Figures  11-4  through  11-6.  A 
highpass  filter  with  a corner  frequency  of  0.8  Hz  was  applied  to  the  accelera- 
tion data  prior  to  evaluation  of  the  spectra.  In  this  case  a spectral  corner 
frequency  of  6.  0 Hz  can  be  seen  and  the  low  frequency  spectral  level  was 
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FIGURE  II-4 

NORMALIZED  ACCELERATION  S-WAVE  ENERGY  DENSITY 
SPECTRUM  OF  N45E  COMPONENT  OF  STATION  2; 
EARTHQUAKE  OF  FEBRUARY  7,  1974 
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FIGURE  II-5 


NORMALIZED  VELOCITY  S-WAVE  ENERGY  DENSITY 
SPECTRUM  OF  N45E  COMPONENT  OF  STATION  2; 
EARTHQUAKE  OF  FEBRUARY  7,  1974 
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FIGURE  II - 6 

NORMALIZED  DISPLACEMENT  S- WAVE  ENERGY  DENSITY 
SPECTRUM  OF  N45E  COMPONENT  OF  STATION  2; 
EARTHQUAKE  OF  FEBRUARY  7,  1974 
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calculated  as  5.44  x 10*  cm  - sec.  Similar  spectra  were  calculated  for  each 
of  the  other  components  recording  the  event. 

The  results  of  this  analysis  and  the  imnlied  source  character- 
istics are  given  in  Table  11-2.  As  would  be  expected  for  an  event  of  this  low 
magnitude,  the  seismic  moments  and  dislocation  radii  are  of  low  values.  How- 
ever, as  this  analysis  does  not  consider  factors  such  as  attenuation,  the  cor- 
ner frequencies  given  are  lower  bounds  and  thus  the  dislocation  radii  could  be 
even  smaller.  However,  attenuation  should  have  only  slight  effect  on  and 

thus  the  estimated  values  of  Mq  are  probably  valid. 

D.  THE  JULY  6,  1974  EARTHQUAKE 

The  July  6,  1974  earthquake  was  recorded  by  six  stations  of  the 
Bear  Valley  accelerometer  array.  In  addition  this  magnitude  3.2  earthquake 
was  also  recorded  in  the  intermediate  field  which  provided  a hypocentral  loca- 
tion and  a fault  plane  solution  for  this  event.  The  epicenter  of  this  event  is 
shown  in  Figure  II- 1 by  the  circle  labeled  E-2  (McEvilly  and  Johnson,  1974). 

A depth  of  6.37  kms  was  determined  for  the  hypocenter.  Both  the  Haskell 
model  and  spectral  estimates  of  the  source  characteristics  were  made  for 
this  event.  The  time  domain  process  will  be  discussed  first. 

The  accelerometer  stations  recording  the  event  were  Stations 
L and  4 which  lie  directly  on  the  surface  trace  of  the  San  Andreas,  3 and  6 
which  lie  to  the  southeast  of  the  fault  and  Sites  7 and  8 which  are  located  on 
the  lower  velocity  northeastern  crustal  block.  The  acceleration  waveforms 
recorded  at  the  various  stations  are  of  a highly  complex  nature  and  it  is  again 
difficult  to  pick  what  might  be  considered  a dominent  pulse  for  use  in  time  do- 
main waveform  fitting  in  the  velocity  traces  of  most  sites.  As  an  example, 
the  acceleration  and  velocity  waves  recorded  at  Stations  3,  Figures  II-7  and 
II-8,  and  7,  Figures  II- 9 and  II- 1 0,  respectively  are  shown.  In  both  cases 
the  high  variability  of  signal  content  frcm  component  to  component  make  selec- 
tion of  a dominant  pulse  difficult. 
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ACCLERATION  (g/IO) 


TIME  (SECONDS) 

FIGURE  II-7 


ACCELEROGRAMS  FROM  STATION  3 FOR  THE 
JULY  6,  1974  EARTHQUAKE 
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TIME  (SECONDS) 

FIGURE  II- 9 

ACCELEROGRAMS  FROM  STATION  7 FOR  THE 
JULY  6,  1974  EARTHQUAKE 
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DERIVED  VELOCITY  TRACES  FROM  STATION  7 
FOR  THE  JULY  6,  1974  EARTHQUAKE 
(DASHED  LINE  REPRESENTS  HASKELL  MODEL  SOLUTION) 
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To  overcome  the  difficulty,  the  whole  space  velocities  calcu- 
lated previously  were  used  to  determine  approximate  P-  and  S-wave  arrival 
times.  The  velocity  of  the  crustal  block  on  which  the  station  was  located 
was  used  for  eacl  station.  For  Stations  2 and  4,  those  lieing  directly  on  the 
fault  trace,  times  based  on  each  model  and  an  average  velocity  were  used. 

As  most  stations  recorded  a strong  P-arrival  on  at  least  one  component,  it 
was  apparent  that  errors  of  as  much  as  one  second  existed  in  these  calcula- 
tions. For  this  reason  the  dominant  pulse  used  for  fitting  with  the  Haskell 
model  was  selected  by  using  the  S-P  times  to  locate  the  approximate  S-arrival 
time  and  then  picking  the  nearest  apparent  dominant  pulse. 

As  previously  explained,  once  the  pulses  to  be  used  in  the  fit- 
ting process  are  selected,  unit  dislocation  waveforms  for  each  station  are 
calculated  using  as  many  external  restrictions  as  possible.  In  this  case  the 
hypocentral  location  and  a fault  plane  solution  have  been  determined  for  the 
source  from  intermediate  field  data  (McEvilly  and  Johnson,  1974).  The  fault 
plane  solution  defined  a fault  striking  N40W  and  nearly  vertical  in  dip.  How- 
ever, this  was  taken  as  the  strike  direction  soley  because  it  parallels  the  San 
Andreas  surface  trace  in  the  vicinity  of  the  epicenter.  The  auxilary  plane  of 
N50E  could  also  define  the  fault  and  thus  both  planes  were  used  to  produce 
theoretical  waveforms  by  Haskell’s  source  model. 

Using  least  squared  error  fitting  on  the  amplitudes  of  the  select- 
ed pulses  from  all  stations  simultaneously  and  separately,  the  initial  source 
structures  were  perturbed  to  produce  7 solutions  defined  by  the  minimum 
residual  over  all  test  cases  for  each  set  of  fittings.  Each  solution  was  thm 
used  to  produce  theoretical  waveforms  at  every  station  which  were  judged  for 
overall  'goodness  of  fit'.  Of  these  solutions  for  the  July  earthquake  the  one 
determined  by  waveform  fitting  at  Site  6 separately  was  selected.  At  best, 
the  quality  of  this  fit  can  be  considered  marginal,  in  general  being  poorer 
than  the  fits  shown  in  Figures  11-8  and  II- 1 0 for  Stations  3 and  7,  respectively. 
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t rom  * igures  II - 8 and  II- 10  it  would  appear  that  the  modeled  pulse  duration 
is  less  than  the  observed  indicating  error  in  the  calculated  fault  length,  rise 
time  or  rupture  velocity.  At  other  stations  the  theoretical  duration  is  closer 
to  the  observed  value.  However,  this  may  indicate  that  tnese  parameters  are 
in  error  to  some  degree.  The  source  parameters  are  given  in  Table  11-3. 

That  the  Haskell  model  could  not  produce  a satisfactory  solution  at  all  sta- 
tions is  probably  the  result  of  several  problems  which  will  be  discussed 
later  in  this  report. 

In  addition  to  the  time  domain  modeling  of  the  earthquake 
source,  spectral  estimates  of  the  source  characteristics  were  also  made  for 
this  event.  As  is  apparent  on  consideration  of  the  empirical  velocity  wave- 
forms in  Figures  II-8  and  II- 1 0 the  spectra  calculated  for  each  station  are 
highly  variable  and  no  single  set  of  spectra  could  be  considered  typical  of 
this  event.  However,  the  S-wave  spectra  of  the  N45E  component  at  Station  8 
are  shown  in  Figures  II - 1 1 to  11-13.  The  corner  frequency,  f^  , for  this  com- 
ponent is  somewhat  lower  than  for  other  stations  but  appears  as  a well  defined 
peak  in  the  velocity  spectra  at  about  4 Hz.  The  results  for  all  stations  are 
shown  in  Table  II-4.  As  with  the  February  earthquake  the  relatively  high  cor- 
ner frequency,  indicating  small  fault  area,  and  low  seismic  moments  calcu- 
lated for  this  event  are  as  expected  from  a low  local  magnitude  earthquake. 

E.  DISCUSSION  OF  RESULTS 

While  the  general  levels  of  seismic  moments  calculated  for 
both  the  February  and  July  earthquakes  are  reasonable  when  compared  to  the 
seismic  moments  and  local  magnitudes  of  other  events  studied  from  the  near- 
zone, as  shown  in  Figure  11-14,  it  is  apparent  that  a high  level  of  variation 
in  spectrally  determined  seismic  moment  exists  between  stations  of  the  near- 
field array  for  each  event.  While  this  is  in  part  accounted  for  by  variations 
in  the  seismic  source  radiation  patterns,  a comparison  of  results  studied  in 


11-21 


i 


T \B  LE  11-3 

BEAR  VALLEY  EARTHQUAKE  OF  JULY  6,  1974 
HASKELL  SOURCE  MO  DE  L SOLU  T ION 


Epicente  r 

Depth 

Velocities 


36  34.  1 'N 
1 2 1 ° 1 0.  1 ' -V 

Southwest  Crustal  Stations: 
Northeast  Crustal  Stations: 


V =5.76  km/  sec 
P 

V = 3.  32  km/  sec 

s 

V =4.56  km/  sec 
P 

V =2.60  km/  sec 
s 


Strike  Angle  : 

N40W 

Dip  Angle  : 

80°  to  NE 

Fault  Dimensions: 

Length  = 0.  18  km 

Width  =0.18  km 

Rise  Time  : 

0.  02  sec 

Rupture  Velocity  : 

2.  5 km/  sec 

Dislocations  : 

Strike-Slip:  6 

Dip-Slip  : 5 

Seismic  Moment  : 

7 n 

8.  3 x 10^  dyne 
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FIGURE  11-11 

NORMALIZED  ACCELERATION  S-WAVE  ENERGY  DENSITY 
SPECTRUM  FOR  N45E  COMPONENT  OF  STATION  8; 
EARTHQUAKE  OF  JULY  b,  1974 
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FIGURE  11-12 

NORMALIZED  VELOCITY  S-WAVE  ENERGY  DENSITY 
SPECTRUM  FOR  N45E  COMPONENT  OF  STATION  8; 
EARTHQUAKE  OF  JULY  6,  1974 
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FIGURE  11-13 


NORMALIZED  DISPLACEMENT  S-WAVE  ENERGY  DENSITY 
SPECTRUM  FOR  N45E  COMPONENT  OF  STATION  8; 
EARTHQUAKE  OF  JULY  6,  1974 
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N45E  4.4  0.62 

N45W  4.0  0.37 

Vertical  6.8  0.06 

Total  Spectra  4.4  0.  72  0.29 
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this  report  and  ot  the  June  22,  1973  Bear  Valley  earthquake  (Turnbull,  et  al.  , 
1979;  Turnbull  et  al.  , 1974b)  indicate  that  the  transmission  and  site  effects 

may  be  more  important  in  causing  these  variations. 

In  Table  11-5  the  seismic  moments  calculated  by  spectral  means 
for  the  three  Hear  Valley  earthquakes  are  shown  for  each  station  recording  an 
event.  It  would  appear  from  this  data  that  the  seismic  moment  found  at  each 
site  is  dependent  on  location  of  the  site.  The  average  moment  for  Stations  7 
and  8,  which  are  located  to  the  northeast  of  the  San  Andreas  Fault,  tend  to 
be  higher  than  for  Stations  2 and  4 which  lie  on  the  fault  and  Stations  3 and  6. 
which  lie  to  the  southwest,  seem  to  have  the  lowest  moments.  McEvilly  and 
Johnson  (1974)  have  found  a similar  correlation  between  seismic  moment  for 
a given  event  and  site  location.  Based  on  their  data  they  find  a factor  of  be- 
tween  3 and  4 difference  between  stations  to  the  southwest  and  to  the  north- 
east  of  the  fault  with  the  stations  on  the  fault  having  a seismic  moment  inter- 
mediate between  these  stations.  In  addition,  however,  it  would  appear  from 
Table  11-5  that  even  for  stations  on  the  same  structure  significant  var.at.ons 
occur.  This  can  be  seen  in  the  consistent  lower  moment  estimate  at  Station 

4 as  compared  to  Station  2. 

Systematic  variations  in  the  acceleration  amplitudes,  as  indi- 
cated by  the  variations  in  seismic  moment,  can  explain  the  difficulty  of  the 
Haskell  model  to  fit  the  July  6,  1974  earthquake.  First,  using  multiple  sta- 
tion fitting,  it  is  apparent  that  one  is  attempting  to  find  the  fault  whose  radta- 
tion  pattern  can  model  the  empirical  amplitude  pattern.  However,  if  the  em- 
pirical data  is  significantly  altered  by  means  external  to  the  source,  as  ap- 
pears to  be  the  case  at  Bear  Valley,  then  it  highly  unlikely  that  a source 
model  can  reproduce  these  effects  with  a 'reasonable'  earthquake  source. 
Secondly,  the  complexities  of  the  structure  produce  anomalous  phases  which 
reduce  the  dominance  of  any  single  pulse  in  the  empirical  wavetratn.  Under 
these  conditions  it  is  highly  unlikely  that  the  same  arrival  will  be  selected  a. 
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TABLE  II - 5 

SEISMIC  MOMENTS  FOR  EACH  EVENT  AT  BEAR  VALLEY 
AS  CALCULATED  AT  EACH  STATION 


6/22/73 

2/7/74 

7/6/74 

Station  2 

1. 8 x 1022 

, 21 

4.  0 x 10 

21 

4.4  x 10 

Station  3 

21 

4.4x10  1 

Station  4 

6.5  x 1021 

4.7  x 1020 

21 

2.1  x 10 

Station  6 

-- 

1.2  x 1 02  1 

Station  7 

22 

5.4  x 10 

-- 

3.  1 x 102‘ 

Station  8 

7.5x10 

3.  2 x 1021 

9 x 1021 

J 


all  stations  used  in  the  fitting  process.  This  in  turn  results  in  incorrect 
phase  and  amplitude  relations  between  the  stations  making  a 'good  fit'  within 
such  restrictions  as  the  fault  plane  solution  nearly  impossible.  Thus,  it  can 
be  concluded  that  the  failure  of  the  Haskell  model  in  this  case  is  not  failure 
of  the  source  model  but  rather  of  its  inability  to  cope  with  transmission  ef- 
fects. 

It  is  also  possible  that  these  transmission  or  site  effects  can 

explain  the  apparent  disagreement  of  seismic  moments  calculated  for  the 

June  22,  1973  earthquake  from  near,  intermediate,  and  far-field  data.  Using 

2 3 

Haskell's  model  a seismic  moment  of  1.2  x 10  dyne*  cm  was  calculated 
for  this  event.  However,  moments  of  much  lower  magnitude  were  found  in 
the  intermediate  and  far-field.  Smith  (1973)  calculated  an  intermediate  zone 
moment  of  2x10  dyne  * cm  while  the  best  far-field  estimates  ranged  be- 

y y 

tween  0.7  and  2 x 1 0 dyne  * cm  (D.  Sun,  personal  communication,  1975). 

If  one  assumes  that  the  variations  of  amplitude  result  from  site  amplifications 

for  stations  on  the  fault  and  to  the  northeast  of  it,  then  using  McEvilly  and 

Johnson's  correction  factor,  the  recorded  amplitudes  of  Stations  7 and  8 should 

be  reduced  by  a factor  of  about  3.  5.  Reducing  these  amplitudes  will  reduce 

the  modeled  dislocations  and  seismic  moment  an  equal  amount.  Thus  a seis- 

22 

mic  moment  of  3.  4 x 10  dyne*  cm  would  be  obtained,  in  much  better  agree- 
ment with  the  intermediate  and  far-field  results  for  this  earthquake. 
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A.  INTRODUCTION 

For  the  past  several  months,  the  examination  of  far-field 
spectra  for  source  characteristics  has  been  conducted  along  several  lines  of 
investigation.  The  investigation  of  two  California  earthquakes  (22  June  1973 
and  28  November  1974)  using  far-field  surface  wave  data  has  been  extended 
from  the  last  study  (Turnbull,  et  al.  , 1975).  In  Subsection  B,  we  discuss 
the  effect  a new  source  region  structure  (McEvilly,  1975)  would  have  on  the 
source  mechanism  solution.  The  investigation  of  theoretical  higher  mode 
spectra  has  also  been  continued,  with  comprehensive  spectral  displays  pro- 
duced and  analyzed  for  two  different  repiesentative  earth  structures  in  Sub- 
section C.  Finally,  in  Subsection  D,  we  present  preliminary  results  on  the 
investigation  of  a swarm  of  earthquakes  that  occurred  in  Sinkiang,  China  in 
1974. 

B.  CONTINUED  EVALUATION  OF  THE  FAR-FIELD  SOURCE  MECH- 
ANISM FOR  TWO  CALIFORNIA  EARTHAUAKES 

In  previous  work  (Turnbull,  et  al.  , 1Q75),  it  was  felt  that  one 
of  the  limitations  to  the  far-field  solutions  for  these  events  was  knowledge  of 
the  source  region  structure.  From  a field  survey  conducted  by  McEvilly 
(1975),  we  were  able  to  obtain  better  structure  estimates  for  both  sides  of 
the  San  Andreas  fault  in  the  Bear  Valley  region.  This  structure  is  listed  in 
Table  III- 1 , with  the  dispersion  curves  for  both  sides  of  the  fault  displayed 
in  Appendix  A.  The  major  difference  between  the  structure  on  each  side  of 
the  fault,  of  course,  is  the  higher  velocity  layers  in  the  structure  on  the 
southwest  side. 
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TABLE  III-  1 

BEAR  VALLEY  REGION  STRUCTURE  (McEVILLY,  1975) 


Depth 

h 

Layer 

Thickness 

Ah 

Northeast  Side  of  Fault 

1 

Southwest  Side  of  Fault  | 

■■ 

V 

s 

P 

V 

P 

V 

s 

0.  6 

0.  6 

2.  0 

2.  38 

1.51 

m 

m 

2.  46 

2.  5 

1.9 

2.  1 

3.  04 

1.76 

2.7 

5.2 

3.  0 

5.  5 

3.  0 

2.7 

5.88 

3.  30 

2.  74 

6.  34 

3.  66 

10.  0 

4.  5 

2.  75 

4.  84 

2.  80 

2.  74 

5.  96 

BB 

15.  5 

5.  5 

2.  75 

4.  98 

2.  88 

2.  8 

6.  94 

4.  0 

19.  3 

3.8 

2.  8 

6.  02 

3.  48 

2.  8 

6.  67 

3.  86 

21.  0 

1.7 

2.  9 

6.  55 

m 

3.  0 

6.  55 

3.  78 

25.  0 

4.  0 

3.  0 

6.  86 

3.  97 

3.  0 

6.  5 5 

3.  78 

50.  0 

25.  0 

3.  32 

8.  2 

4.  65 

3.  32 

8.2 

4.  65 

Merge  With  Gutenbcrg-Bullcn  Earth  Model 
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Another  way  to  determine  the  impact  the  structural  differences 
have  on  the  observed  spectra  is  to  generate  theoretical  spectra  using  a source 
imbedded  in  these  structures.  Using  the  bodywave  source  mechanisms  for 
these  two  events,  spectral  variation  with  respect  to  depth  and  radiation  pat- 
terns of  a particular  depth  were  generated.  The  two  source  mechanisms 
are  listed  below  in  Table  I1I-2. 


TABLE  III  - 2 

SOURCE  PARAMETERS  FOR  THE  TWO  CALIFORNIA  EVENTS 


^ Event  Description 

Depth  - h 
(km) 

Dip  - 8 
(deg) 

Slip  - X 
(deg) 

Strike  - 0 
N (deg)  E 

Bear  Valley  6-22-73 

10.25 

70 

30 

135 

Central  California  11-28-74 

4 

90 

0 

35 

These  spectra  and  radiation  patterns  are  displayed  for  both 
events  in  Appendix  B.  The  major  conclusions  that  can  be  drawn  from  com- 
parisons between  the  two  structural  responses  are: 

• From  the  displays  of  Rayleigh  and  Love  wave  amplitude  spec- 
tra as  a function  of  depth  we  see  that,  for  the  appropriate 
source  depth,  and  for  both  events,  the  spectral  levels  for  the 
northeast  structure  are  approximately  an  order  of  magnitude 
higher  than  those  levels  generated  by  the  southwest  side. 

• From  the  Rayleigh  and  Love  wave  radiation  patterns  for  both 
events  generated  at  20  seconds  period,  we  see  that  the  pat- 
terns have  shapes  that  are  structure  dependent  as  well  as 
having  different  spectral  levels. 
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With  these  structural  properties  firmly  in  mind,  we  attempt- 
ed to  'fit'  the  observed  spectra  in  the  same  manner  as  in  previous  analysis 
of  these  events  (Turnbull,  et  al.  , 1975).  The  stations  available  for  the  Bear 
Valley  event  of  22  June  1973  are  displayed  in  Figure  III-  1 , with  the  Rayleigh 
wave  traces  used  from  OGD  and  ALQ,  and  the  Love  wave  trace  from  LASA. 

In  the  spectral  fitting  process,  no  absolute  minimum  residual  was  observed, 
but  a solution  was  obtained  by  picking  the  local  minimum  residual.  Because 
the  stations  lie  to  the  east  of  the  fault,  the  northeast  structure  was  used. 

The  solution  was  found  to  be  close  tc  the  bodywave  mechanism,  with  a depth 

of  6 km,  dip  of  80°,  slip  of  0°,  strike  of  N175°E,  and  a seismic  moment  of 
22 

0.  9 x 10  dyne-cm.  This  moment  is  about  an  order  of  magnitude  less  than 
that  obtained  by  fitting  near-field  acceleration  data,  and  approximately  the 
same  as  that  obtained  previously  using  the  Helmberger  structure. 

Analysis  of  the  central  California  event  of  28  November  1974 
presented  a more  complex  situation.  The  stations  available  are  displayed  in 
Figure  111-2,  with  the  Rayleigh  wave  traces  used  from  NORSAR,  TLO,  ALP  A, 
ZLP,  and  KIP,  and  the  Love  wave  trace  from  LASA.  With  stations  on  both 
sides  of  the  fault,  a problem  arises  as  to  which  source  region  structure  is 
appropriate.  An  exhaustive  approach  was  applied  to  this  situation,  with  solu- 
tions obtained  by  using  the  northeast  model  only,  the  southwest  model  only, 
and  for  the  first  time,  a combinacion  of  the  two  structures.  The  solutions 
obtained  for  the  different  combinations  'if  stations  and  structures  based  on 
the  minimum  residual  criteria  and  the  distribution  of  minimum  residual  cri- 
teria are  given  in  Tables  111-3  and  II1-4,  respectively.  All  of  the  solutions 
yield  a similar  source  mechanism;  a vertical  strike -slip  fault  with  a strike 
direction  of  N45  E.  This  agrees  very  well  with  the  bodywave  solution.  How- 
ever, except  for  the  solution  obtained  by  fitting  the  three  stations  in  the  west 
with  the  southwest  model,  all  of  the  focal  depth  estimates  were  too  deep. 
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FIGURE  111-2 

TRAVEL  PATHS  TO  THE  AVAILABLE  STATIONS  FOR  THE 
NOVEMBER  28.  1974  EARTHQUAKE 
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From  the  above  analysis,  we  have  reached  the  following  con- 
clusions : 

• Fitting  the  surface  wave  spectra  with  one  earth  model  using 
stations  on  both  sides  of  the  fault  neglects  the  effect  of  the 
medium  discontinuity.  One  possible  way  to  correct  for  this 
discontinuity  would  be  to  make  the  appropriate  travel  path 
correction.  This  problem  will  be  examined  in  future  studies. 

• Fitting  the  surface  wave  spectra  with  two  earth  models  would 
seem  to  be  a reasonable  approach,  if  we  consider  the  effect 
of  the  discontinuity  as  a difference  in  the  source  excitation. 
One  problem  in  this  approach  is,  however,  that  the  combina- 
tion of  the  two  earth  models  for  simulating  both  sides  of  the 
fault  actually  neglects  the  possible  physical  complexity  of  the 
source  due  to  the  discontinuity. 

• At  the  present  time,  we  have  no  solution  to  the  above  pro- 
blems. Since  the  location  of  this  event  is  definitely  on  the 
west  side  of  the  San  Andrea:  fault,  the  far -field  solution  ob- 
tained  by  using  the  three  stations  to  the  west  of  the  fault  with 
the  southwest  earth  model  should  be  considered  a good  esti- 
mate. 

C.  CONTINUED  EXAMINATION  OF  THE  THEORETICAL  HIGHER 

MODE  SPECTRA 

Continuing  with  our  examination  of  theoretical  first  higher 
mode  Rayleigh  and  Love  wave  spectra,  double  couple  sources  were  placed 
at  three  representative  depths  (4,  30,  and  50  km)  in  a Gutenberg-Bullen 
earth  model  and  a Hamilton-Healy  earth  model.  Variations  cf  the  spectra 
with  respect  to  dip,  slip,  and  st-ike  were  examined  (see  Appendix  C).  We 
can  summarize  the  properties  of  this  spectra  as  follows: 
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• Regardless  of  the  structure,  the  only  consistent  'holes'  in  the 
spectra  occur  for  Rayleigh  waves  generated  by  a vertical 
strike-slip  source  mechanism.  In  particular,  for  the 
Gutenbe r g - Bullen  structure,  two  'holes'  appeared  for  this 
source  mechanism  placed  30  km  deep  at  9 and  20  seconds 
period. 

• The  level  of  first  higher  mode  spectra  is  approximately  an 
order  of  magnitude  less  than  that  of  the  fundamental  mode 
spectra  over  the  same  frequency  range  (15  to  40  seconds  per- 
iod). The  higher  mode  spectra  over  the  3 to  15  second  period 
range  is,  however,  approximately  the  same  level  as  the  fun- 
damental mode  15  to  40  second  spectra. 

• For  the  Gutenberg-Bullen  structure,  the  first  higher  mode 
spectra  drops  off  rapidly  for  periods  greater  than  9 seconds, 
especially  for  Love  waves.  This  drop-off  is  even  more  severe 
for  the  Hamilton-Healy  structure,  starting  at  6 seconds  per- 
iod. This  is  a direct  result  of  the  lower  velocity  material  in 
the  upper  crust  of  the  Hamilton-Healy  model.  Therefore,  we 
should  expect  to  observe  higher  modes  with  greater  regular- 
ity from  thick  high  velocity  source  region  crustal  structures 
using  presently  installed  instrumentation. 

D.  ANALYSIS  OF  THE.  S1NK1ANG,  CHINA  EARTHQUAKE  SWARM  OF 
AUGUST  1974 

During  most  of  the  month  of  August  and  the  early  part  of  Sep- 
tember 1974,  a major  earthquake  swarm  of  over  300  recorded  events  occur- 
red in  Sinkiang,  China.  The  location  of  the  swarm  is  shown  by  the  boxed 
area  in  Figure  111-3  (This  map  was  part  of  an  earlier  study  by  Turnbull, 
et  al.  , 1974a).  The  number  and  magnitude  of  ^hese  events  as  a function  of 
date  is  shown  in  Figure  111-4. 
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FIGURE  111-3 

LOCATION  OF  S1NK1ANG,  CHINA  SWARM  OF  AUGUST.  1974  - (BOXED  AREA) 
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FIGURE  III -4 

THE  NUMBER  AND  MAGNITUDE  CHRONOLOGY  OF  THE 
RECORDED  EVENTS  OF  THE  SINKIANG  EARTHQUAKE 
SWARM  OF  AUGUST  1974 


This  swarm  was  chosen  for  intensive  analysis  because  of  the 


I 
I 

following  reasons: 

• The  opportunity  to  analyze  a large  number  of  events  from  a 
small  region  - investigate  the  repeatability  of  source  mech- 
anism from  event  to  event  and  local  tectonics. 

• The  possibility  of  using  a master  evert  or  events  with  spec- 
tral ratios  to  measure  path  effects. 

• The  opportunity  to  determine  effective  attenuation  by  eliminat- 
ing  radiation  pattern  effects  (knowing  the  source  mechanism) 
and  by  two  station  path  corrections  (in  the  instances  where 
they  occur). 

Preliminary  analysis  of  five  of  these  events  have  consistently 
yielded  shallow  source  depths  in  the  4 to  10  km  range.  The  source  mech- 
anism details,  though,  have  not  followed  a particular  trend  at  this  point  of 
the  study.  A more  detailed  analysis  of  this  swarm  will  be  presented  in  the 
next  report. 
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SECTION  IV 

THE  M - m RELATIONSHIP 
s b 


In  previous  work  by  Turnbull  et  ai.  (1974b,  1975)  the  Mg-  rr»b 
discriminant  was  investigated  with  the  objective  of  reducing  the  scatter  of  the 
earthquake  population.  Several  methods  were  tried  to  achieve  such  a reduc- 
tion; all  of  them  designed  to  modify  in  various  ways  the  determinations  at 

individual  single  stations  of  the  VLPE  network,  and  assess  the  impact  on  net- 
work magnitude  measurements.  However,  no  signilicant  reduction  of  the 

M - m scatter  was  observed, 
s b 

In  this  report  we  focus  upon  the  determination,  and  apply 
the  maximum  likelihood  method  of  Ringdal  (1975)  in  an  attempt  to  improve  the 
accuracy  of  the  POE  bodywave  magnitudes  for  our  data  base.  The  basic  idea 
is  that  PDE  magnitudes  of  small  events  tend  to  be  biased  high,  because  only 
a few  stations  (with  particularly  favorable  radiation  patterns)  report  such 
events,  and  network  magnitudes  are  therefore  obtained  by  averaging  over  a 
set  of  stations  that  is  not  representative  of  the  whole  network.  By  taking  into 
account  the  estimated  noise  levels  at  the  non-detecting  stations,  this  magni- 
tude bias  can  partly  be  compensated  for.  Our  model  is  as  follows: 

• For  large  events,  the  typical  number  of  PDF  stations  that  re- 
port a bodywave  magnitude  is  assumed  to  be  at  least  Nq  (where 

N is  a number  to  be  estimated) 
o 

• For  a given  event,  the  number  N of  stations  that  actually  re- 
port a magnitude  is  compared  to  N_.  If  N > Nq,  it  is  assumed 
that  essentially  all  PDE  stations  detected  the  event  If  N < Nq, 
we  assume  that,  for  each  of  (N.-N)  stations,  the  magnitude  was 
less  than  a certain  threshold  value  m . 


1V-1 


• The  threshold  values  m at  each  station  are  assumed  to  be  in- 

2 

dependency  sampled  from  a normal  distribution  m~N(m  , rj  ) 

o o 

• The  world-wide  magnitude  values  of  the  event  are  assumed  to 

2 

follow  a normal  distribution  N(/i,(7  ) ; where  the  oarameter  M 
denotes  the  true  magnitude  of  the  event. 


The  likelihood  function  (Ringdal,  1 97 5 ) becomes;  for 


Mm)  = 


where  the  m.  denote  the  magnitude  measured  at  the  detecting  stations,  and 
where  and  0 are  the  standard  normal  density  and  distribution  functions, 
respectively.  An  estimate  of  n is  then  obtained  by  numerical  maximization 
of  (IV-  1 ). 

In  order  to  do  this,  it  is  necessary  to  have  available  realistic 

estimates  of  the  parameters  <7_,  m , a , N ; besides  the  observed  quantities 

S o o o 

N and  | m.;  i = 1 N }.  Studies  by  Evernden  (1969)  of  the  noise  levels  of 

some  of  the  better  WWSSN  stations  seem  to  indicate  values  of  m = 5.  0 and 

o 

a = 0.  2 as  realistic  estimates.  Veith  and  Clawson  (1972)  found  a = 0.4  to 
o S 

be  representative  of  the  WWSSN  network.  From  studies  of  PDE  bulletins,  we 

have  found  a value  of  N = 25  as  a reasonable  estimate  of  the  number  of  sta- 

o 

tions  expected  to  report  magnitudes  of  large  Eurasian  events. 

Figure  IV-1  shows  the  theoretical  bias  (m  - M)  between  the 
average  magnitude  m (over  the  N detecting  stations)  and  the  true  magni- 
tude [I  as  a function  of  the  ratio  N/N  and  for  various  values  of  (m  - m ). 

o o 

Equation  (IV-1)  has  been  used  to  compute  these  bias  values.  Note  that  the 

only  assumed  parameter  values  in  Figure  IV-1  are  <7  = 0.4  and  a = 0.  2; 

S o 


N < N : 
o 


IV  -2 


J 


Ratio  N / N 


FIGURE  IV-1 

EXPECTED  BIAS  OF  CONVENTIONAL  MAGNITUDE  m RELATIVE 
TO  TRUE  MAGNITUDE  [l  FOR  A HYPOTHETICAL  NETWORK  AS 
A FUNCTION  OF  THE  PERCENTAGE  OF  DETECTING  STATtqnS 


hence  the  graphs  can  also  be  applied  to  networks  with  different  parameter 

values  N and  m than  used  he  re. 
o o 

Table  1V-1  lists  the  events  in  the  data  set  along  with  certain 
event  parameters.  Among  those  are  the  NORSAR  and  PDE  reported  rr>b 
values;  the  number  N of  PDE  station  magnitudes  available,  and  the  maxi- 
mum likelihood  PDE  magnitudes  from  the  proceeding  model  Figures  IV-2 
and  IV -3  plot  NORSAR  magnitudes  versus  PDE  standard  m^  and  PDE  maxi- 
mum likelihood  (M-L)  m^,  respectively.  It  is  seen  that  the  latter  case  pro- 
duces a linear  trend  with  a slope  close  to  1.  00,  while  Figure  IV-2  exhibits 
the  network  magnitude  bias  effect  for  small  events,  and  produces  a more 
horizontal  slope  of  about  0.  8.  The  conclusion  is  that  the  maximum  likelihood 

m values  are  more  consistent  estimates  of  Ji  than  the  standard  PDE  mfa  for 
b 

our  data  sot. 

In  order  to  see  if  the  use  of  maximum  likelihood  ir>b  values 

would  reduce  the  scatter  in  the  Mr-  mb  population,  V J_PE  was  plotted 

against  PDE  nv  (Figure  1V-4)  and  against  PDE(M-L)  m (Figure  1V-5).  The 
b 

straight  line  of  each  figure  is  a least  squares  orthogonal  fit,  and  the  orthogon- 
al standard  deviation  Cr  has  been  computed  together  with  the  slope  « and  in- 
tercept (1  of  the  estimated  functional  relationship: 


M = a • m.  + P 
s b 


(IV  -2) 


We  note  that  the  value  of  a does  not  decrease  when  PDE  m^s  are  replaced 
by  maximum  likelihood  estimates;  on  the  contrary,  a slight  increase  is  found 
(from  0.  238  to  0.262).  However,  the  slope  a is  much  steeper  in  the  first 
case  (1.34  versus  1.23),  hence  the  values  of  O are  not  directly  compatible. 
(For  example,  an  estimation  method  that  assigned  identical  values  to  all 
events  would  produce  a vertical  slope  with  CT  = 0;  hence,  a low  value  of  n 
does  not  necessarily  imply  improved  estimation  techniques.) 
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In  order  to  examine  more  closely  the  variatiors  of  the  M - m^ 
slope  as  a function  of  estimation  techniques  of  M ^ and  m^,  a total  of  nine 
cases  were  run  based  on  the  given  event  population.  In  each  case,  the 
values  estimated  by  either  V LPE  (averaging),  ALiPA  or  NORSAR  were  com- 
bined with  the  m values  of  NORSAR,  PDF,  and  PDE(M-L).  The  results  are 
b 

summarized  in  Table  1V-Z.  The  following  points  are  noteworthy 

• Four  of  the  runs  produce  a virtually  identical  slope  («  ~ 1.  37). 

These  are  precisely  those  four  runs  that  combine  'consistent' 

M and  m values,  i.  e.  , values  free  of  network  bias  effects, 
s b 

(NORSAR  or  ALPA  versus  NORSAR  or  PDE(M-L)  mfa.  ) 

• A consistently  high  value  of  the  slope  (1.66  or  1.64)  is  found 

when  PDE  m is  plotted  against  a consistent  M . 

b s 

• Aconsistently  low  value  of  « (1,23  or  1.24)  results  when  V LPE 
M is  plotted  against  a consistent  m^. 

• When  PDE  m,  is  plotted  against  V LPE  M , n is  again  high, 

b s 

showing  that  the  network  bias  effects  in  PDE  magnitudes  dom- 
inate those  of  VLPE. 


I 
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l 

I 
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Hence,  the  behavior  of  the  computed  slope  agrees  well  with 
what  could  be  expected  from  network  bias  consideration,  and  it  appears  that 
the  most  accurate  linear  functional  relationship  between  and  m^  for  the 
given  data  set  (ranging  in  values  from  about  4.  0 to  6.  0)  has  a slope  of 
approximately  1.4. 

Considering  closer  the  four  cases  of  consistent  slope  estimates, 

it  is  interesting  to  note  that  the  value  of  cr  is  lower  when  using  PDE(M-L)  in 

versus  either  ALPA  or  NORSAR  M compared  to  when  NORSAR  m is  used 

s o 

(a  x 0.  26  versus  a ~ 0.  31).  It  would  be  interesting  to  compare  PDE(M-L) 
to  the  VLPE- ALPA- NORSAR  combined  network  with  the  network  bias  reduced 
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detection  at  these  stations. 


by  maximum  likelihood  processing;  however,  we  have  not  been  able  o do  " is, 
mostly  because  ol  .he  lack  ot  reliable  V LPE  data  lor  some  stations. 
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APPENDIX  A 

DISPERSIVE  CHARACTERISTICS  OF  THE 
BEAR  VALLEY  STRUCTURE 


Using  the  velocity  profile  obtained  by  McEvilly  (1975),  dis- 
persion curves  were  generated  for  the  structures  representing  both  sides  of 
the  fault.  Figure  A-l  shows  the  velocity-density  profile  for  the  northeast 
side,  and  the  associated  Rayleigh  and  Love  wave  dispersion  curves  are  dis- 
played in  Figure  A-2.  Similarly,  the  velocity-density  profile  and  associated 
dispersion  curves  are  given  in  Figures  A-3  and  A-4,  respectively. 

The  major  difference  between  these  structures  is  the  higher 
velocity  material  on  the  southwest  side  of  the  fault.  This  is  reflected  in  the 
dispersion  curves,  with  higher  group  and  phase  velocities  between  10  and  30 
seconds  period  for  the  southwest  structure. 


h (km) 


I 


Same  as  Gutenberg-Bullen  Earth  Model 


FIGURE  A- 1 

BEAR  VALLEY  STRUCTURE  - NORTHEAST 
SIDE  OF  FAULT 
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Velocity  (km/ sec) 


APPENDIX  B 

SPECTRAL  RESPONSE  OF  THE 
BEAR  VALLEY  STRUCTURE 


Theoretical  spectra  were  generated  for  each  event  using  the 
Bear  Valley  structures  as  describ  -d  in  Appendix  A.  Of  particular  interest 
were  changer  in  spectral  shape  and  level.  For  the  Bear  Valley  event  of  22 
June  1973,  Figures  B-la  and  B-lb  give  the  Rayleigh  and  Love  wave  ampli- 
tude spectra  as  a function  of  depth  for  the  northeast  structure,  while  Fig- 
ures B -2  a and  B-2b  give  similar  spectra  for  the  southwest  structure.  Fig- 
ires  B-3a  and  B-3b  show  the  Rayleigh  and  Love  wave  radiation  patterns  at 
20  and  30  seconds  period  for  the  northeast  structure,  and  Figures  B-4a  and 
B-4b  give  the  patterns  for  the  southwest  structure. 

For  the  central  California  event  of  28  November  1974,  Fig- 
ures B-5a  and  B-5b  give  the  Rayleigh  and  Lowe  wave  amplitude  spectra  as 
a function  of  depth  for  the  northeast  structure,  while  Figures  B-6a  and  B-6b 
give  similar  spectra  for  the  southwest  structure.  Figures  B-7a  and  B-7b 
show  the  Rayleigh  and  Love  wave  radiation  patterns  at  20  and  30  seconds 
period  for  the  northeast  structure,  and  Figures  B-3a  and  B-8b  give  the  pat- 
terns for  the  southwest  structure. 
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SOURCE  MECHANISM  OF  THE  22  JUNE  1973  BEAR  VALLEY 
EVENT  IN  THE  NORTHEAST  STRUCTURE 
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FUNDAMENTAL  MODE  LOVE  WAVE  SPECTRA  USING  THE  SOURCE 
MECHANISM  OF  THE  22  JUNE  1973  BEAR  VALLEY  EVENT 
IN  THE  NORTHEAST  STRUCTURE 
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FUNDAMENTAL  MODE  RAYLEIGH  WAVE  SPECTRA  USING  THE 
SOURCE  MECHANISM  OF  THE  22  JUNE  1973  BEAR  VALLEY 
EVENT  IN  THE  SOUTHWEST  STRUCTURE 
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FUNDAMENTAL  MODE  LOVE  WAVE  SPECTRA  USING  THE  SOURCE 
MECHANISM  OF  THE  22  JUNE  1973  BEAR  VALLEY  EVENT 
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APPENDIX  C 

THEORETICAL  FIRST  HIGHER  MODE 
SURFACE  WAVE  SPECTRA 


As  in  previous  work  (Turnbull,  et  al.  , 1975),  theoretical  first 
higher  mode  Rayleigh  and  Love  wave  spectra  were  generated  using  the  Hark- 
rider  (1964)  formulation  of  a double  couple  source  in  a layered  half  space. 
Our  objective  is  to  identify  any  special  spectral  characteristics  (such  as 
shape  or  amplitude  level)  that  would  be  useful  in  determining  source  depth 
or  structure  in  the  source  region. 

Figures  C-la  through  C-ld  display  the  variation  of  dip,  slip, 
and  strike  for  first  higher  mode  Rayleigh  and  Love  waves  over  the  period 
range  of  3 to  15  seconds  and  for  three  representative  source  depths  in  a 
Gutenberg-Bullen  earth  model.  ULing  the  Hamilton-Healy  earth  model,  the 
effect  of  dip  angle  variation  is  displayed  over  the  same  period  range  and  for 
three  source  depths  in  Figure  C-2.  Finally,  in  Figure  C-3,  higher  mode 
surface  wave  spectra  over  the  period  range  of  15  to  40  seconds  with  dip 
angle  variation  is  shown  for  three  source  depths  in  a Gutenberg-Bullen 
earth  model. 
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ABSTRACT  (ConTiimr  on  revrrtt  side  II  ntcetfry  rind  Identity  by  block  numbrr)  “ — — 

Several  continuing  investigations  of  the  seismic  source  using  near- 
and  far-field  data  are  discussed.  We  completed  our  analysis  of  near-field 
acceleration  data  recorded  at  Bear  Valley,  California,  by  examining  two 
small  events  which  occurred  on  7 February  1974  and  6 July  1974.  The 
February  event,  using  only  spectral  analysis,  was  found  to  have  a moment 
of  3 x 10  dyne-cm  with  an  equivalent  circular  dislocation  radius  of  0.  18 
km  . Both  spectral  and  time  domain  waveform  fitting  techniques  were 
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20.  continued 


applied  to  the  July  event.  The  scdj>Hons  using  each  method  were  in  rea- 
sonable agreement,  with  a seismic  moment  in  the  range  of  2 to  4 x 1 02  ^ 
dyne-cm  over  an  equivalent  dislocation  radius  of  about  0.  1 to  0.  3 km2. 

Using  far  - field  surface  wave  data,  two  central  California  area 
earthquakes  (22  June  1973  and  28  November  1974)  were  re-examined  using 
two  new  earth  models  which  were  determined  by  McEvilly  (1975)  for  the 
northeast  and  southwest  side  of  the  San  Andreas  fault.  Several  far-field 
solutions  for  the  source  mechanism  were  obtained  by  using  various  com- 
binations of  stations  and  earth  models,  including  use  of  both  structures  at 
once  for  the  first  time.  However,  for  the  central  California  event  of  28 
November  1974,  the  mechanism  which  agreed  closely  with  that  obtained 
by  other  means  used  only  those  stations  on  the  west  side  of  the  fault, 
with  the  southwest  structure  used  in  the  spectral  fit. 

Theoretical  first  higher  mode  surface  wave  spectra  were  generated 
for  a double  couple  source  in  a layered  half  space.  Both  Gutenberg -Bullen 
and  Hamilton-Healy  earth  models  were  used.  As  with  fundamental  mode 
spectra,  'holes'  were  found  almost  exclusively  for  vertical  strike-slip 
faults.  Higher  mode  spectral  levels  in  the  3 to  9 second  period  range  were 
found  to  be  comparable  with  15  to  50  second  period  fundamental  mode  spec- 
tra. Also,  from  the  excitation  of  the  earth  models,  greater  amplitude 
higher  mode  signals  from  shallow  events  should  occur  from  high  velocity 
crustal  source  regions. 

Preliminary  results  have  been  obtained  for  five  large  events  which 
were  a part  of  »he  Sinkiang,  China  earthquake  swarm  of  August  1974.  Both 
simple  spectial  and  spectral  ratio  fitting  methods  were  employed.  The 
far-field  soiitions  showed  a consistent  trend  for  these  five  events  of 
shallow  focai  depth. 

Finally,  we  concluded  our  examination  of  the  scatter  of  the  earth- 
quake population  in  the  Mg-  rn^  discriminant.  Focusing  on  the  bodywavc 
magnitude  estimate  instead  of  the  surface  wave  magnitude  as  in  previous 
studies,  we  applied  the  maximum  likelihood  method  of  Ringdal  (1975)  in 
an  attempt  to  improve  the  accuracy  of  the  PDE  bodywave  magnitudes.  Al- 
though the  scatter  of  the  population  dia  not  change  significantly,  the  slope 
of  the  data  set  decreased  approximately  0.25. 


UNCLASSIFIED 


